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APPLICATION OF MAXIMUM LIKELIHOOD TECHNIQUES TO 
THE DESIGN OF OPTIMUM STAR TRACKERS 


By Edwin C. Foudriat 
Langley Research Center 

SUMMARY 

An analytical method has been developed for the design of optimum star trackers. 

In the method both the measurement procedure - that is, the position where measure- 
ment samples are made - and the measurement processing - that is, the estimation of 
the star parameters from the measurement samples — are optimized. The method for 
measurement processing is based upon the use of maximum likelihood techniques. The 
equations to be satisfied for estimation of the position and intensity of a star are given 
with the assumption that the photoelectrons leaving the cathode obey Poisson statistics. 

With the properties of the maximum likelihood method and the Cramer -Rao 
inequality, a technique is developed for optimization of the measurement procedure. 
General equations are derived for the optimum measurement position on the cathode. 

The procedure is applied to the design of a slit -type star tracker where it is 
desired to determine the x position of a star image with a Gaussian intensity distribution. 
The results indicate the star -tracker slit width should be based upon star image size and 
be independent of star-to-background intensity ratio. The measurement procedure indi- 
cates sample should be made along the limb of the star image. In addition, the measure- 
ment processing equation for the single-slit star tracker can be simplified so that an 
approximate estimate of the star position offset from its assumed position can be obtained 
in closed form. Results indicate the solution error to be small and dependent upon star- 
to-background intensity ratio. 


INTRODUCTION 

In the determination of the attitude of spacecraft, stellar observations have been 
found to be both extremely useful and accurate. Many spacecraft, especially those like 
the ones discussed in references 1 and 2 which require precision attitude determination, 
have used star-tracker techniques. Many methods of observation and many mechaniza- 
tions of the star-sighting technique have been employed. The accuracy of the techniques 
has been predicted analytically and verified experimentally. However, to date no work 


was found which indicates how measurements might be taken in an optimum fashion and 
how these methods might be used to design an optimum star tracker. Therefore, an 
effort was made to apply optimization techniques to the design and mechanization of a 
star-tracker system. 

As has been well established, the arrival of stellar energy can be described by a 
statistical process. Therefore, the analysis immediately employs the statistical proce- 
dure called point estimation (ref. 3) to predict the position of the center-of-stellar 
energy. The technique of point estimation most frequently used is maximum likelihood, 
and it is by this method that a minimum variance estimator is derived. This leads to 
an optimum method of processing the measured information. The analysis is carried 
still further to minimize the variance function which in turn results in an optimum 
design. To illustrate the principle, the procedure is carried out for the specific 
stellar energy distribution represented by an error function which should be generally 
representative of the intensity of starlight when focused through a well-designed optical 
telescope. Hence, the result of this analysis is both the optimization of the star-tracker 
parameters and the procedure for processing the measurement to obtain the estimate of 
star position. 


AA 

B 

C 

E (> 

E n 

F b 

F s 


h 


f 2 


SYMBOLS 


area on the photocathode sampled by the measurement 


constant, fj(s) f 2 (k) 
A c 

constant, 


x c 


expectation operation 
expectation of number of electrons 

background distribution function integrated over the area AA 
star distribution function integrated over area A A 
function of s defined by equation (31) 
normalization function defined in equation (30) 
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function defined by equation (32) 
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normalized point background distribution function 
normalized point star distribution function 
constant specifying the background intensity 
nondimensional slit width in x-direction 
likelihood function 

mean intensity per unit area 

number of samples 

number of photoelectrons 

number of photoelectrons in ith sample 

general statistical density function 

square of x relative position from center-of-star intensity, (x 
value of q where optimum measurement should be taken 
function defined by equation (24) 
function defined by equation (26) 

slit width in y-direction 
time interval of measurement 
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time interval of ith measurement 
measurement position on photocathode 
position of center of star image 
assumed position of image prior to measurement 
maximum photocathode dimensions 
position of ith measurement 

ratio of photoelectron rate of ith measurement to expected background 
rate 

constant, 4 s — 

Kb 

displacement of true image position from its assumed position 
y relative position from center -of-star intensity, y - y 

parameter 

value of 6 which maximizes L 

intensity of background, photoelectrons/unit time 

intensity of star, photoelectrons/unit time 

function of random samples which define parameter estimate 

value of (J. which maximizes L 

ratio of expected signal rate to expected background rate 


4 



empirical constant 


r dummy variable 

APPLICATION OF MAXIMUM LIKELIHOOD PROCEDURES TO 
STAR POSITION AND MAGNITUDE DETERMINATION 

It is assumed that the optical system of a telescope focuses its gathered light upon 
a photocathode placed at the focal plane of the system. As had been indicated by many 
authors, this photon collection is a statistical process so that the number of photoelec- 
trons leaving the cathode at the point x,y can be given in terms of its probability dis- 
tribution function. It is well established that, for sufficiently long integration times, the 
photoelectrons leaving the cathode obey Poisson statistics (ref. 4). 

A pictorial representation of the star image and background intensity 
(photoelectrons/unit time) distribution is shown in figure 1. The photocathode surface 
is represented by the plane with dimensions x c and y c . The magnitude of the inten- 
sity is shown by the amplitude normal to the x,y plane. A star positioned at x,y with 
intensity distribution A s G s (x,y,x,y) and a background distribution given by A^G^CXjy) 
is illustrated. The quantities A s and A b represent the total photoelectron emission 
rate resulting from the star and background, respectively, where the factors affecting 
collection efficiency, effective aperture size, and photocathode efficiency have all been 
taken into account in the term A. The functions G s (x,y,x,y) and G b (x,y) are the 
normalized intensity density functions so that 

G(x,y) dx dy = 1 (1) 

Before proceeding with the analysis it should be noted that the pictorial represen- 
tation of figure 1 shows the expected situation - that is, if a number of unit time samples 
of intensity were made at the point x,y, the mean of these samples would be 

m = A s G s (x,y,x,y) + A b G b (x,y) (2) 

and that the samples would obey a Poisson distribution given by 

p (n) = m^ (n = 0,1,. . .) (3) 

In addition, it is assumed that the samples are statistically independent. 
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The samples discussed so far are measurements at a specific point. In general, 
measurements must be taken over an area of the photocathode and over a specific time 
interval. Thus, each sample is given by 


A s F s (x,y,x,y)At = X s \ \ G g (x,y,x,y) dA At 

A b F b (x,y)At = A b f f G^y) dA At 
J J AA 


(4) 


where AA specifies the area sampled at the photocathode and At the time interval of 
the measurement. 

If measurements of the photoelectron emissions are taken at various points on the 
photocathode, then with equation (3) the joint probability density for N samples can be 
written as 


N 

p(n 1; . . .,n N ) = J7 p(ni) 
i=l 

N 


- JT •jjA s F s (xi ,Yi,x,y) + ^b F b( x i’yi)] 


ni exp i,yi,x,y) + X b F b (x i ,y i )]At i 


n i ; 


(5) 


Equation (5) can be rewritten as 


/ v \ n |[As F s( x i>yi> x >y) + A b F b( x i>yi)] At i} ni 

)(n^,. . . ,n^ = exp(- ^ |j^ s F s (x b y b x,y) + ^- b F b (x b yjJ| AtjJ j]^ — 1 ” 


H- 


This probability density function describes mathematically the probability that measure- 
ments made of the photocathode current at the N points x b y^ will result in photo- 
electrons at the ith measurement for all N measurements. 

A statement of the problem to be solved is as follows: when the n^ measurements 
are given, estimate the position of the central point of the star image and the star inten- 
sity - that is, what values should be assigned to x,y and A s . This general class of 
problem is given in reference 3. A procedure employed for the solution of this type of 
problem is maximum likelihood. The procedure can be described by considering a ran- 
dom sample nj,. . .,njq- from a distribution having a density function given by p(n,0) 
where 6 is the parameter to be found. Then, the likelihood function of the random 
sample can be defined as 
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It is obvious that the likelihood function is identical to the joint probability density func- 
tion. The procedure to maximize L is to determine a function ju(nj,. . .,n^ such 
that when 


e = p(n 1} . . .,n N ) (8) 

L is a maximum. The maximum statistic - that is, the value of 6 which actually 
maximizes L for the measurements n^,. . .,n^ — is designated by 

§ = p(nj,. . .,n N ) ( 9 ) 

The maximization of L is accomplished by the standard procedure of partial differenti- 
ation with respect to the parameters of interest and subsequent solution of the resultant 
equations . 

Applying this procedure to the problem at hand, equation (6) can be used with equa- 
tion (7) to give 


L(A.s,x,y,n].,. • -,n N ) = p(nj,. . .,n N ) 

Taking the log of both sides gives 

L’^sA^nj,. . .,n N ) = - ^ [x s F s ( Xi , yi ,x,y) + X b F b ( Xi , yi ) At t 

i 

+ ^ n i 1 ogJ|x s F s (x i ,y i ,x,y) 
i ^ 

+ ^Fb^yJlAtij - ^ log ni ! 


(10) 


( 11 ) 


Taking the partial derivatives of L' with respect to the parameters of interest gives 


— = - y f s Ati + y ni 7 — — — T = o 

8X S i i ( X sF S + ^b F b) 


— = - y x s au ^ + y ^ 

9x A s 1 3x Lj 1 


^s 


a_Fg 

9x 


3x Y (A S F S + ^b F b) 


= 0 


(12) 


(13) 



(14) 


9L’ 

9y 


l 


Xr At,- 


3F S 

9 y 


2- 


0 Fr 

Xs ~8y~ 
( x s F s + A b F b) 


= 0 


The simultaneous solution of equations (12) to (14) for A g , x, and y provides the solu- 
tion for the extreme conditions of L'. Determining A s , x, and y which maximize 
L' gives the maximum likelihood function estimation of the star magnitude and loca- 
tion, respectively. A more general discussion of the above problem can be found in 
reference 5. 

A few general assumptions can be made concerning the representation of the inten- 
sity rate distribution function of the signal and background and the general measurement 
procedure. First, it is assumed that the background rate distribution is constant; that is, 

*b F b( x >y) = X x b y A = K b ( 15 ) 

c^c 


This assumption can be made without loss of generality since all the intensity rate varia- 
tion at the photocathode can be adequately described in the function F s (x,y,x,y). Second, 
the assumption that At is equal for all measurements is reasonable. Employing these 
assumptions, equation (12) can be written as 

= 0 < 16 > 

Kif 



The factor n^/At is the photoelectron arrival rate. Thus, 

n i 1 _ _ Photoelectron rate of measurement 

At 1 Expected background arrival rate 


and 


where 


A s F s + j _ Expected photoelectron rate _ j + p f x y £ 
Kb Expected background rate \ ’ ’ ’ 


. . Expected signal rate 

p(x,y) = “ - - 

Expected background rate 

Heuristically, equation (16) indicates that the value X s is selected which makes the 
measured and expected arrival rates match, and the function F s is used as a means of 
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weighting the various measurement points. If it is assumed that a single measurement 
is made, then 


K b(“l ~ !) 
s F s 


(17) 


under the assumption that F g (x,y) ^ 0, and x,y are known exactly. 

Equation (13) can be processed in a similar manner. Using these factors in equa- 
tion (13) gives 



(18) 


The only difference between equations (18) and (16) is that the weighting function on the 
individual terms is altered from F s to dF s /dx. 

In general, it is obvious that a trial-and-error procedure can be used to obtain the 
solution to the parameters A g , x, and y with equation (16), equation (18), and an 
equation similar to equation (18) for y. Equations (13) and (14) have identical form; 
therefore, the result is applicable to both. 


ESTIMATION OF VARIANCE AND METHODS FOR THE DETERMINATION 
OF THE OPTIMUM MEASUREMENT PROCEDURE 


The maximum likelihood procedure has certain statistical properties which permit 
estimates of the accuracy with which the parameter 8 is determined. These properties 
are known as efficient and sufficient estimates. These can be stated as 

1. If an efficient estimate of the parameter 9 exists, the likelihood equation will 
have a unique solution for 8. 

2. If a sufficient estimate of the parameter 9 exists, any solution of the likelihood 
equation will be a function of 9. 

A detailed discussion of sufficient and efficient estimates is given in reference 6. One 
of the properties of estimations is that the variance satisfies the Cramer-Rao inequality 
which for an unbiased estimator is given by 


Var(e) § 



(19) 


For the case of an efficient estimator, the equality sign can be used (ref. 6). 
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An additional property of the maximum likelihood procedure is described as follows 
Under certain general conditions, the likelihood equation has a solution which converges 
in probability to the true value of 6 as n — °°. This solution is an asymptotically nor- 
mal and asymptotically efficient estimate of 9. Hence, if enough sample measurements 
are taken, the equality sign in equation (19) can be used. 

It is apparent from equation (19) that minimization of the variance is related to 
maximizing the absolute value of E n {d%L ' The procedure for maximizing 
E n { 92 L'/ 9 0 2 } can be carried out by 


9E r 

9x 


= 0 


9E 


n 


9 y 


= 0 


( 20 ) 


Using x as an example and starting with equation (13) and equal measurement times 
gives 


9 2 L* 

9x 2 




5U^X B At 
+ p dk 2 s 


_ £ 9 F S ^ Xs At 


i 8x (1 + P) 2 


(21) 


From the definition of a it is apparent that 

E n{«i} = 1 + P 

and equation (21) upon taking the expectation becomes 


( 22 ) 



AtA s 2 

^9F s \ 2 

K b 

V 9x / 


1 + P 



i 


(23) 


It is apparent that |E n {92L'/ 9x 2 }| 
maximized. Mathematically, if 


is maximized if each term of the series is 


Rxfoy) 


AtA^/9FA 2 
E b \ 8 x / 

1 + p 


(24) 


then 


£Rx 

9x 


= 0 


9R- 


x 


9 y 


= 0 


(25) 
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can be used to fulfill the criteria of equation (20). The value of x,y obtained from this 
solution can be used to determine the position on the focal plane where the measurements 
should be taken in order that the optimum-*- prediction of the parameter x can be made. 

A similar procedure can be carried out for y giving an equation identical to equa- 
tion (24) and for A s giving 


AcAtFs 2 
%s = K b (l + p) 


(26) 


where the maximization is carried out by use of equation (26). 


EXAMPLE OF STAR-SIGHTING PROCEDURE 
FOR SINGLE PARAMETER CASE 


The procedures outlined in the previous two sections can be illustrated by an 
example using a single parameter. The example is simple enough to carry through mathe- 
matically yet has some practical utility so that the results can form a basis for suggesting 
the procedure in more complex cases. 

For the example, it is assumed that G s (x,y) is an error function intensity distri- 
bution at the focal plane given by 


G(x,y,x,y)= i exp(- 


(x - x) 2 (y - y) 2 
2 + ~ 2 ' 


(27) 


The error function curve used to represent the star image has been normalized. This 
means that a unit distance (x - x) on the focal plane is related to the condition where the 
star intensity has dropped to 0.606 of its maximum value. Hence, the analysis is based 
upon the normalized distance along the focal plane which is related in turn to the optical 
image diffraction pattern (blur circle), the f-number, and aperture of the star tracker. It 
is assumed in the example that any telescope distortion is reasonably well corrected 
along both x- and y-axes so that the diffraction pattern is both circular and symmetrical. 

A comment concerning the validity of the error function representation is in order. 

It is well known that for monochromatic plane wave light incident upon a circular aper- 
ture, the diffraction pattern of an optically perfect telescope is a circular Bessel function. 
In this pattern, distinct rings (ref. 7) are present with the distance of each ring having a 
definite relationship to the wavelength of the source. When the source is no longer mono- 
chromatic but consists of a wide spectrum, then the rings are no longer distinct, either 


^Optimum in the sense of minimum variance. 
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becoming blurred or losing the character of rings altogether. Because stellar detection 
of position uses a wide spectrum of light, it is felt that the error function represents a 
realistic approximation to the intensity distribution at the focal plane and hence a good 
choice for the analysis. 

The first step is to determine the signal intensity from a finite aperture at the 
focal plane. Assume a reactangular aperture 2s -units wide in the y-direction and 
k-units wide in the x-direction. Equation (4) can be written as 

F s ( x ,x,s,k) = i exp ^ * exp(- (28) 

2 

Equation (28) assumes that y is known and that the slit is positioned such that 
£ = y - y. Hence, the £ integral is symmetrical. By using the substitution that 
rj = t - x, equation (28) can be rewritten as 

x+ — 

F s (x,x,s,k) = J j* exp(-^)d^ ^ ^ exp (-^ dr] (29) 

0 \ « _ y _ " 


The error integral is tabulated and published in numerous places, and therefore the inte- 
gral of equation (29) is easily obtained numerically. A plot of the normalized function of 
the second integral is shown in figure 2 for values of k ranging from 0.2 to 4.0. Fig- 
ure 2 indicates all solutions have the same shape so that an analytical representation of 
F s will be sought in the form 

F s (x,x,s,k) = fj (s) f 2 (k) f 3 (x - x,a x ) (30) 

From equation (29), it is obvious that 

fl<s) = » ip exp ('T") d5 <31> 


For fixed slit width in the y-direction, fj(s) can be readily obtained from a set of nor- 
mal probability distribution tables. The function f 3 (x - x,cr x ) is estimated to be of the 
form 


f 3 (x-x,a x ) = exp 


f (X - x) 2 ' 

V , 


(32) 


In figure 2 is a curve of equation (32) for cr x = 2.20, and this curve shows almost perfect 
agreement with the curve of the normalized function F g (x,x,s,k) for k = 1. In order to 
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obtain other values, cr x (k) was calculated at x - x = 1.0 for various values of k. A 
plot of this function is shown in figure 3. The final function f 2 (k) is the normalization 
factor; that is, the value of 


f "' X+ 2 / 

L , * exp v 


o' 

_ 222 


x-x-l 


dr; 


when x - x = 0. The function f 2 (k) j <j2ii is plotted in figure 4. 

With equations (30) and (32) it is a relatively simple matter to obtain explicit values 
for the functions used in the previous section. Thus, p(x,x) and 9F S //9x become, 
respectively, 


p(x,x) 


^s fj(s) f2(k) 


K 


b 



(x - x) 2 \ 
°x ) 


(33) 


and 


9x 


= B exp| 


(X - x)‘ 




2(x - x) 


J x 


(34) 


where B = fj(s) f 2 (k). The maximum likelihood solution (eq. (18)) becomes 


I 


-1 + 




X S B 

K b 


f 3 (x - x) 


(x - x) f 3 (x - x) = 0 


(35) 


As indicated in the analysis of the previous section, the heuristic approach to the 
solution of equation (35) is to make each term in the summation equal to zero (or as close 
as possible) while the weighting function on each is the term (x - x) fg(x - x). The 
terms where x = x and f 3 (x - x) ~ 0 result in a zero weighting which means any mea- 
surement either at peak signal intensity or where the signal intensity is zero contains 
very little information about x and hence are of no value. 

Equation (24) can be used to construct the optimum measurement. technique. Sub- 
stituting equations (33) and (34) gives 


Rx - 


4 AtX s 2 B 2 (x - x) 2 (-2(x - x) 2 

Kb 0 


°x 2 


“'X 


1 + B exp 
Kb 


(x - x) 2 


(36) 
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Using equation (25) and algebraic manipulation gives 


4 AtX s 2 B 2 

9R x _ K b°x 2 
ax 


2(x - x) f 3 2 + ^(x - x) f 3 3 + 2(x - x) 2 f 3 f 3x + ^(x - x) 2 f 3 2 f 3? 


1 + l^ f3 ^ X " ^ 


(37) 


where 


f 3 = 

a x a v 


(38) 


Setting the term in the numerator of equation (37) equal to zero and using equation (38) 
give 


(x - x) f 3 2 


1 + 


A s B ( x - x) 2 / X S B 


K b 


*3 - 


2 + 


J x 


K b 


= 0 


(39) 


The previous discussion indicated that no information is available when x = x and 
f 3 (x - x) = 0, since the weighting function at these positions was zero. Hence, x = x and 
f 3 2 = 0 must be terms which minimize R^. The term within the brackets in equa- 
tion (39) satisfies the maximum condition and its solution is given by 


(x - x)2 - + P(X ’^ 

2 + p(x,x) 


(40) 


If q = (x - x) 2 , equation (40) can be written by using equation (33) as 


where 


q = 


1 + /3 exp(-^ 


2 + )8 exp(-i 


/3 = 


A s B 

K b 


(41) 


(42) 


The value of q which satisfies equation (41) gives the position at which the mea- 
surement should be taken in order that the variance of the estimate will be minimized. 
This solution will be designated q. Once q is known, equations (19) and (23) can be 
used to calculate the minimum variance as 


14 



Var(x) = 


J x 


2 + ^ ex p(-^ 


4 AtNK b j3 2 exp (- ^ 


(43) 


where N is the total number of measurements made at the point q. 

Review of equation (30) demonstrates that the solution of q and Var(x) are 
dependent upon k and s, the x- and y-dimensions of the slit, respectively. Therefore, 
the optimum measurement position and the minimum variance depend upon the slit con- 
figuration. Hence, the solution to the design of the instrument is to seek that value of 
(k,s) which reduces equation (43) to its minimum value. 

For the example considered, the optimization of the slit width is carried out in the 
x-direction only by assuming that the slit in the y-direction extends the full length of the 
photocathode; that is, s = y c » 1. Hence, equation (42) becomes 


= x s f 2(k) _ c f2(k) 
" k 

x c 


(44) 


where K b is given by equation (15) and the value of fi(s) (eq. (31)) is 
equations (41) and (43) can be rewritten as 


— Lr. Thus, 
^7T 


Ox(k) 


q = 


f,(k) / Q 

1 + C 4 — exp ( — 


°x( k )y 


2 + cM^exp 


e x ( k ) 


and 


Var(x) = 


'jM 

.NAtxJ 


e x (k) 






4C 


f 2 2 (k) 


exp 


2q 

CT x(k), 


(45) 


(46) 


Equation (45) was solved numerically to obtain q as a function of slit size for 
various values of C. The values are plotted in figure 5 where \jl[, the optimum mea- 
surement position, is plotted against k, the slit width. Since C is the ratio of star 
intensity to background intensity per unit distance, then figure 5 shows the effects of 
relative changes in the star-to-background intensity ratio. 
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Examination of figure 5 indicates that for very small slit widths the optimum 
measurement position lies between 1.0 and 1.4. As the slit width increases, the posi- 
tion ql/2 increases. For example, when C = 100 and k = 2.0, the edges of the slit 
fall at 0.62 and 2.62. It is readily apparent from examination of figure 2, which can be 
interpreted as normalized star intensity as a function of slit position for various slit 
widths, that optimum measurements are made near the maximum slope of the star 
intensity. This is logical from the information theory standpoint since the slit is posi- 
tioned to avoid the high intensity at the center of the image which provides considerable 
increase in noise with little increase in information as to the star position. As the star- 
to-background intensity ratio decreases, the optimum position moves inward. The 
inward movement is a result of trading background signal, which results in increased 
noise only, for increased star image intensity, which results in a slight increase in 
information along with the noise increase due to higher intensity. For very low star 
intensities, the inward shift ceases. 

Using the solution of equation (45) for q, it is possible to calculate the minimum 
variance for each case. Figure 6 shows the plot of equation (46) for the function 

y ar (x) as a function of slit width for various values of the parameter C. Since 

\l2ir 

the Var(x) is a direct function of C, the curves have been plotted with a scale factor 
for convenience. It is clear that for optimum variance, independent of intensity ratio, a 
slit width of about 2.0 gives the lowest value on each curve. Thus, for the optimum 
rectangular-slit star tracker, which has an error function image intensity shape, its slit 
should be twice its image size. Image size is defined as the point where the amplitude 
G s is equal to 0.606 of its peak value. 

PROCEDURE FOR DETERMINING STAR POSITION BASED ON 
OPTIMUM MEASUREMENT TECHNIQUES 

The preceding section has determined the optimum star -tracker design. It is of 
interest to determine how one might use the maximum likelihood formulation to obtain an 
optimum estimate of the star position. The starting point of such a technique is equa- 
tion (18), which must be satisfied for the measurements made. The procedure in the 
previous section (eq. (45)), indicates that two points symmetrically spaced on each side 
of the maximum star intensity point x can be used if known. 

The procedure for determination of the maximum intensity point can be illustrated 
in the sketch shown as figure 7. The position of the image assumed prior to the measure- 
ment is shown as the solid curve with x a as its center. The center of the slit is placed 
alternately at xj and x 2 which are ±(q distance from x a > and the time increment 


16 



at each measurement position is equal. The true image position, shown as a dashed line 
in figure 7, with center at x is displaced a distance e from the assumed position. 

The procedure to be developed is a determination of the value of e using equation (18) 
as a starting point. 

From equations (33) and (34) 


_ As tlW m exp (- fa - it 

K b V °x 


and 


8F C 


3x 


r (xi,x) = B 
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-*x 


exp - 
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v^2N 
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where i = 1,2. By using 


x 2 = x a + r 
X 1 = x a " r ) 


(47) 


where r = ]fq, equation (18) can be written after algebraic manipulation as 
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where /3 is given by equation (44) . 
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is difficult. Therefore, 
e « r, and r ~ then 

(49) 


With these approximations, equation (48) can be solved for e as 


(50) 


Equation (50) has been checked computationally against the solution of equation (18) 
for the condition of equal measurement on each side of the image circle. The results of 
these calculations are shown in figure 8. Comparison of the true e with that calculated 
from equation (50) for values ranging from 0.001 to 1.0 indicates that the error is a fixed 
percentage of the true value and depends upon the factor C. Thus, the approximation 
appears to be correctable by an empirically derived function dependent upon star inten- 
sity, and hence a useful means of obtaining control error information for the star tracker. 

CONCLUSIONS 

From the analysis given in this report, the following conclusions can be made: 

1. Assuming convergence to an efficient estimate of variance, the Cramer-Rao 
equation can be used as a procedure to determine the optimum star -tracker design 
parameters. For the example used in this paper, the optimum slit width i£ twice the 
image intensity parameter size and relatively independent of the star-to-background 
intensity ratio. 

2. In addition, the procedure can be used to determine the measurement position on 
the star image in order to obtain the minimum variance estimate of star position. A 
measurement procedure for the optimum slit size which places the center of the slit 
aperture at about 1.5 units from the center position of the image for stars which are 
bright relative to their background is the optimum measurement position. This value 
decreases slightly as the background intensity relative to star intensity increases. 
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Obviously, the solution of equation (48) explicitly for e 
solution by approximation is attempted. If it is assumed that 


expf-2*^tl?\ * expf-^ 


cosh * 1 


I 


3. The maximum likelihood procedure for an optimum two-position measurement 
scheme can be approximated to obtain a simple equation for the error in star position 
from its assumed value. The solution error is small and, in addition, is a fixed per- 
centage of the solution which is dependent upon the star -to -background intensity ratio. 
Hence, further correction by a factor dependent upon star intensity should be possible. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 13, 1967, 

125-19-03-12-23. 
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Figure 7.- Sketch of star-measurement technique. 



Estimated star image offset 



Figure 8.- Comparison between true star position and star position estimated by equation (50). 
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